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Abstract 

In this paper we study Lifshitz tails for continuous Laplacian in a 
continuous site percolation situation. By this we mean that we delete a 
random set from M. d and consider the Dirichlet or Neumann Lapla- 
cian on D = M. d \ F UJ . We prove that the integrated density of states 
exhibits Lifshitz behavior at the bottom of the spectrum when we 
consider Dirichlet boundary conditions, while when we consider Neu- 
mann boundary conditions , it is bounded from below by a van Hove 
behavior. The Lifshitz tails are proven independently of the percola- 
tion probability, whereas for the van Hove case we need some assump- 
tion on the volume of the sets taken out as well as on the percolation 
probability. 
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1 Introduction 



The present study deals with the behavior of the integrated density of states 
for Laplacians in a site percolation setting on the continuum. 

We start by defining the main object of our study, the integrated density of 
states (IDS). The IDS is a concept of fundamental importance in condensed 
matter physics. Indeed, it is considered to measure the number of energy 
levels per unit volume, below a given energy. For P, a finite-dimensional 
orthogonal projection, we denote by tr(P), the dimension of its range. Let 
P(-oo,e] De the spectral projection of a random Schrodinger operator H w and 
Ai be a cube of side length L around the origin. The restriction of P(_oo,e] 
to the cube Al has tr(xA L P(-oc,E]) as a dimension range. Here for a set A 
we denote by xa, the characteristic function of the set A. We denote by \E\, 
the Lebesgue measure of E. We consider 

N{E) = lim -±—tr( X A L P(-oo, E]). (1.1) 

It is called the integrated density of states of H w ( IDS). See [10] for 
alternative definition of the IDS. 

The question we are interested in here concerns the behavior of N at the 
bottom of the spectrum of H u . 

In 1964, Lifshitz [15] argued that, for a Schrodinger operator of the form 
H u = —A + V u , there exists Ci,c 2 > such that N(E) satisfies the asymp- 
totic: 

N(E) ~ ci exp(— c 2 (-E' — Eq)~$), as E \ E . (1.2) 

Here E is the bottom of the spectrum of H u . The behavior ( II. 2p is known as 
Lifshitz tails. In the last thirty years, there has been vast literature, both 
physical and mathematical, concerning Lifshitz tails and related phenomena. 
We do not try to give an exhaustive account of this literature. The work of 
Kirsch and Metzger [TT] gives a survey of such results and basic references 
on this subject. Below, we give results on the IDS behavior in the context of 
percolation. 
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There is a long history of works which consider Hamiltonians on percolation 
graphs. These graphs are obtained by removing sites (site percolation) or 
bonds (bond percolation) from a graph, for example Z d , in a random way. 
For example, one might remove sites in the graph Z d independently of each 
other with a probability p which is the same for each site, in other words there 
is a site at a given i G Z d with probability q — 1 — p. The percolation graph 
(consisting of the sites not removed) has various connected components, the 
so called clusters. It is known that there is a critical value p c which depends 
on the dimension d, such that for q < p c there are only finite clusters, while 
for q > p c (the percolation regime) there is also an infinite cluster. For d > 2 
the critical value satisfies < p c < 1. The above results are discussed and 
proved in [5], they were mainly obtained by Hammersley in the late fifties. 
We remark that the fact, that there is only one infinite cluster, was proven 
thirty years later by Aizenman, Kesten and Newmann, see [S]. 

Quantum Hamiltonian with a percolation structure were first studied by 
de Gennes et al in [3j 0] , where the Hamiltonian of a binary alloy was con- 
sidered. It is proved that the spectrum of these percolation Hamiltonian is 
pure point if the fraction q (of conducting sites) is less than the critical value 
p c . For the percolation regime it is argued in [3] that the spectrum contains 
a continuous part. 

In [2] the authors considered the site percolation model on the lattice 
Z d , for d > 2. They investigated the density of states for the tight-binding 
Hamiltonian projected on the infinite cluster. It is shown that, almost surely, 
the IDS is discontinuous on a set of energies which is dense in the band. See 
also |20j for this type of questions. 

There is also a collection of papers concerning the asymptotic behavior 
of the density of states for Laplacians on percolation graphs (see e. g . pE],[6], 
13 , [12] , [IS] and [H]). For example, the papers [T2] and [TF] investigate the 
asymptotic behavior of the integrated density of states Nd(E) and Nn(E) 
near the bottom E of the spectrum of the bond percolation Laplacian on 
Z d with Dirichlet repectively Neumann boundary conditions. These authors 
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prove that Nu(E) has Lifshitz tails with Lifshitz exponent | for any < 
q < 1. The Neumann Laplacian on the percolation graph has an infinitely 
degenerate eigenvalue at the bottom E = of its spectrum and its density 
of states is discontinuous there. So, it is reasonable to consider N(E) = 
Nn(E) — iVjv(O). The results from p2] and [T7] show that for q < p c the 
quantity N(E) has Lifshitz behavior with Lifshitz exponent |, independent 
of the dimension. For q > p c N(E) has a van Hove behavior near E = 0, i.e. 
N{E) ~ CEt 

In this paper, we consider quantum percolation problems on M. d . In one 
particular model we remove a set Si = S + i near the point i 6 Z rf with prob- 
ability p. More precisely, let {Ci}iez d be a collection of independent {0, 1}- 
valued random variable with P(£ = 1) = p. We set D = M. d \ U&=i ^ an d 
denote by and iff? the Laplacian restricted to D with Neumann and 
Dirichlet boundary conditions respectively. As the set D is random, the op- 
erators Hp and H£ are random operators. 

We note that the set D may be a connected set for all u if the set S is 
small, e.g. if S C A a = {x G lR d | —a < Xi < a for z = 1 . . . d} with a < 1. On 
the other hand, if S is big and p is small, the set D is a union of bounded 
connected components. 

We consider the behavior of the integrated density of states for these 
two classes of operators. We will obtain Lifshitz tails for the Dirichlet case 
under mild assumptions. For the Neumann case we prove that N(E) > CEi 
if the set S is not too big. The result on the Neumann case is based on a 
deterministic result (Theorem 14. ip which we believe is of some interest on 
its own. We apply it to other percolation models as well, for example to a 
Poisson percolation model. 

2 Model and results 

Let us first define the model we are studying. 
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2.1 The Model 



Consider the probability space Q = {0, 1} Z , which is endowed with the 
usual product sigma-algebra, generated by finite dimensional cylinder sets, 
and equipped with a product probability measure P. Elementary events in 
Q are sequences co = (w 7 ) 7gZ d, we assume their entries to be independently 
and identically distributed according to a Bernoulli law. i.e 

F{u = 1}= P . (2.3) 

Here p is the site probability, a parameter in ]0, 1[. 

Let Z d be the lattice, considered as a subset of R d . We set Al the hypercube 
defined by A = [— k, ^[ d - Let S be a bounded subset of R. 
We set 

r«= U ^(7 + 5) = |J (7 + 5) 

where the notation 5 = and 1 S = S was used. 

Let D u = M d \r LJ , be the domain obtained from M. d by taking off the 
random set r w . We denote by the restriction of the Laplacian —A to 
D w , with Dirichlet boundary condition and by the restriction of —A 
to with Neumann boundary condition. and are self-adjoint lin- 
ear operators on L 2 (IR d ), we call them Dirichlet respectively Neumann 
Laplacian for site-percolation in M. d . We denote by No, respectively 
by N N the IDS of H° respectively of H%. 

Remark 2.1. As in JM/, we can write formally, in the Anderson form 
with some random potential i.e 

H° = -A + V u , 

Kj = 0J 1 f{x — 7). with f = 00, on S and elsewhere . 

Let us consider the map V, from fl, to the set of the self-adjoint operator 
on L 2 {R d ), such that to u associate V(co) = H u = H* = — A'^ d ^ r ^,» G 
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{N,D}. V is measurable. Sometimes we suppress the superscript D and 
respectively N and write only H w . In this cases the statement remains true 
for both M D and H" 

Let Ti be the unitary translation operator on L 2 (Z d ), i.e 

%ip{x) = i){x -i), G L 2 {R d ) and x G M. d . 

As the probability measure P is ergodic with respect to the group of trans- 
lation (Ti) i£ i d i acting as %{oj) = (w 7+i ) 7eZ d, we get 



By this, we deduce that H u is a measurable family of self-adjoint and ergodic 
operators. According to [El El HE] we know that there exists E, E pp , E ac and 
E sc closed and non-random sets of R such that E is the spectrum of if^ 
with probability one and such that if a pp (respectively a ac and a sc ) denote 
the pure point spectrum (respectively the absolutely continuous and singular 
continuous spectrum) of H u , then E pp = a pp , E ac = a ac and E sc = a sc with 
probability one. 

The following Lemma gives the precise location of the spectrum. 
Lemma 2.2. The spectrum E, of is [0,+oo[ with probability one. 
Proof: First let us notice that for any oj G f2, we have 



(2.4) 



H u >0. 



(2.5) 




[0, +oo[c E for P — almost every a> G fi. 



(2.6) 



For this, let f2, be the following events 




oj G f2 : for any fceN, 



There exists A£ C IT, such that D 



R' 
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Here A.( u ) is the set of points which are both in A and A£. In (12. 7\\ we 

k 

asked that all sites inside of to be present. Let E G [0,+oo[= E(— A) 
be arbitrarily fixed. Using Weyl criterion, we know that there exists a Weyl 
sequence {^E,n)nen C L 2 {R d ) , for —A. Thus ||^_E, n ,|| = 1, for all n G N and 

lim \\{A + E-I)cp E>n \\ = (2.8) 

Notice that for any i G Z d , {Ti^E^n&i is also a Weyl sequence. Without loss 
of generality, we assume that the sequence (<PE,n)neN is compactly supported. 
So for any uj G Q, there exists a Weyl sequences {(fiE,n)nen for (—A) on R d 
with the property that all the supports are contained inside the cubes of 
(12.71) . So for every u G f2 and any n G N there exists an integer and a 
cube A^J and as in (12. 7p such that supp(ip% n ) is contained in A ^; . That 
is, 

min{| x-y\:xe supp^ n ; y G R*\Ajj$} > 0. 
So, for any n G N and u G Q, we get 

ll(^-^)^JII = ll(A + i?-I)^,J|. (2.9) 

Hence, (vlrJneN is also a Weyl sequence for f^. So we get ( 12. 6ft for any 
a; G fi. Now it suffices to check that P(f2) = 1. For this let A be an integer 
bigger than 2. (A fc A ) AeN C M d be a sequence of disjoint cubes in R d i.e 
A Ml n A kM = whenever A x ^ A 2 . We set fl fcjA = {u G : DJ£ M = R^J . 
So (fifc,A)AeN, is a sequence of two by two statistically independent sets, with 
non-zero probability and independent of A G N. So, Using the Borel-Cantelli 
lemma, we get that P(Ofc) = 1 for any fcGN + 2, where 

P(fik) = limsupfi fciA . 

A— >+oo 

The proof of Lemma 12.21 is ended by noting that 

□ 
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2.2 Main results 

Our study is on the bottom of the almost sure spectrum of H u . We recall 
that | S | denotes the volume of S 

Theorem 2.3. For p e]0, 1[, we have 

e->o+ \ogE 2 v 1 

Remark 2.4. We notice that the set D w may have unbounded connected 
component depending on the shape of the set S and on the probability p. 
Indeed the result of the above theorem is independent on the existence or 
nonexistence of such unbounded clusters. 

Theorem 2.5. If \S\ ■ p 1 < — is a bounded set in M. d , we get 

e^o+ log(.E) - 2 v ' 

Remark 2.6. • In fact we prove a slightly better result then $2.11\) . In- 
deed we prove that 

N{E)>C-E%, (2.12) 



with C a constant which can be computed from our proof below. 

In the generality of Theorem \2.5[ we can't bound N(E) from above, as 
we discuss at the end of the paper. 

The result given in Theorem \2.5\ is deduced from a more general result 



of Theorem 4-1 which gives a deterministic lower bound of the numbers 



of eigenvalues for Neumann operators. 



3 The proof of Theorem 12.3 

In this section, we prove Theorem 12.31 
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3.1 Preliminary 

We start by recalling the following result and giving some properties of the 
IDS. 

Proposition 3.1. Let ip G C^°(R d ), then 

lim J-trMHJxAr) = E(tr(xA 1 ^)XA 1 )), (3.13) 

for ^-almost all u. Here E ; is the expectation with respect to the probability 
measure P. 

Proof: First we write xa l — J2ieA L nz d XAi(i)> Here Ai(i), is the cube of 
center i end side length 1. We set d = tr^tp^H^XA^i))- So Q is an ergodic 
sequence (with respect to Z d ) of random variables. So 

L. t r (^(#J XA J = _L J2 &■ (3-14) 



i<EA L r\Z d 



By the Birkhoff's ergodic theorem, the sum in (I3.14p converges to its expec- 
tation value. This ends the proof of ( 13.13!) . □ 
Now, we notice that both sides of (13.131) . are positive linear functionals on 
the bounded, continuous functions. So, they define positives measures and 
H respectively, i.e 

/ ip(\)dfi L (\) = J-tr(<p(H u )xA L ) 

JR \ A L\ 

and 

(p(\)dn{\) = E(tr(xA x ^(iL)XA x ; 



For those two measures we have the following result proved in [20] . 
Theorem 3.2. For any ip G C^°(1R) and for almost all u G £1 we have 

lim (ip,d/j, L ) = (<p,d/j). 

k— >-oo 
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Remark 3.3. We call the non-random probability measure jx the density of 
states measure. It satisfies the following fundamental properties 



N(E) = fi((-oc,E}), 
E(H U ) = supp(jj). 

Let HP (oj) be the operator H® (u) restricted to with Dirichlet bound- 
ary condition also on 8Al , while (u) when we consider Neumann bound- 
ary condition on DAl- As H u is an elliptic operator, H\ (oj), • G {D, N}, has 
a compact resolvent . So H^ (oj) has a discrete spectrum with possibility of 
accumulation at the infinity [T9]. Let us denote the sequences of eigenvalues 
by 

El(A L ) < E'(A L ) < ■ ■ ■< E'(A L ) < ■ ■■ (3.15) 

For E G R we denote by N '(H^ (oj) , E) the number of eigenvalues of H\ l (oj) 
less or equal to E, of course counted with their multiplicities. 
To prove Theorem 12.31 we prove a lower and an upper bounds on N(E). The 
upper and lower bounds are proven separately and based on the following 
result (Theorem 5.25 p. 110 of [18J). 

-±-E{N(H° l (oj),E)} < N(E) < * E{JV(J^>), E)}. (3.16) 
I Al I I Al I 

Inequalities in (I3.16P are based on Akcoglu-Krengel subergodic theorem. In- 
deed 

<M©<<< UA >) (3.17) 

and 

#£ua>) < R°M (3-18) 
hold on L 2 (A 1 U A 2 ), for all bounded cubes A l5 A 2 C M d whenever A 1 n A 2 = 

iia- 



io 



3.2 The Dirichlet case: 
3.2.1 The upper bound 

The upper bound is proved by a comparison procedure. Indeed, let, 

Vu=^2 UjXs(x ~ 7), 



and 

H u = -A + V u 

For any cl d , we set 

Qk L M) = (<P,H<fll>), cp,ijeH 1 (A L \T ul )=V AL , (3.19) 

and 

Qt L M) = (<P,H U 1>), <P,*l>e ^o(Al) = V Kl . (3.20) 
From [19], we recall the following result, 

Lemma 3.4. For any L,k G N*, we have 

sup inf _ Qk L (M) < 



tp u ... ,(p k _ 1 ev AL 1>£[<pi>- ,Vfc-i]- L nD A 



= 1 



sup inf Qt L {w)- (3-21) 

From Lemma 13.4^ one deduces that for any 11 G N*, we have 

E^H^Al)) < E n (H u (A L )). (3.22) 

Thus, we get that for any Bel, 

N(H Al (u),E) < N(H Al (cj),E) (3.23) 

We notice that for it is already known that it exhibits Lifshitz tails, by 
the result of Kirsch and Simon [13] . This ends the proof of the upper bound. 

□ 
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3.2.2 The lower bound 

Let iff be the operator with oj 7 = 1 for any 7 G Z d . We recall that for 
any C l d , we have 

H u < H° l {u) < H° Kl . (3.24) 
So by the min-max argument we get that 

E?(A L ) < E?(H 1M ). (3.25) 
Using equation ( 13.161) one gets, 

N(E) > ±- d .F{E°(A L )<E} 

> — -. ■ F{E? (A z ) < E and V7 G A L n Z d , w 7 = 0} (3.26) 

> — ■ P{£?f (# ,aJ < # and V7 G A L n Z d , w 7 = 0}. (3.27) 
If L is such that 

E?{H q , Kl ) < E, (3.28) 

then 

(ElZD > • P{u, = 0}l A ^l = (1 - V f d . (3.29) 

We remark that (I3.28P is satisfied for L « £?~^. This ends the proof of the 
lower bound. 



4 The proof of Theorem 12.5 



4.1 The deterministic results 

The aim of this section is to prove the following theorem 
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Theorem 4.1. Let Q be a subset ofM. d . Let H be the self adjoint operator 

defined as the Laplacian operator, —A restricted to L 2 (M> d \Q) with Neumann 

boundary conditions on dQ. Assume there exists a sequence of cubes Al such 

2, d 

that for L big enough we have —7 | Ql |< 1 (Ql = AlC\Q). Then, then there 

L d 

is a constant C, such that for E > small enough the IDS N of H N satisfies 

lo g N N (E) A 
e^o+ log£ _ 2 ( J 

Proof: For the proof of we start by recalling the following Lemma from 
1141: 



Lemma 4.2. Let n G N and ipi, ■ ■ ■ , <p n G 7i be in the domain T> of a self- 
adjoint operator A, which acts on a separable Hilbert space H. Suppose there 
are constants a% < • • • < a n < a such that 

\(<fi, <fj) - 5ij\ < ex and \((p h Aipj) - ajdij < e 2 (4.31) 

for all i,j — 1, • • • If ei < 1, then 

N(A,^^)>n. (4.32) 
1 — £1 

Before giving the proof of Lemma 14.21 let us use it to end the proof of 
Theorem 14.11 Let us denote by H\®, the operator H N restricted to A^, 
with Dirichlet boundary condition on 8Al. H^® is the self adjoint operator 
associated with the following quadratic form 

U% D (u,v)= I Vu(x) • Vu{x)dx- (4.33) 
Ja l \u l 

with domain 

V{UD = if e H\A L \n L ); f\dA L = 0}. 

7i^ D is a densely-defined quadratic form. We denote by H^® the self-adjoint 
operator associated to "Hff D . 
For n G N d , we set 



d 



«n(^) = (-) 2 ^n 



2 
i=l 
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and 

^A^) = (j) d/2 fl^), (4.34) 
i=i 

where 

tpk = cos(k7rx); k = 1, 3, 5, • • • , . 
We note that the family (3> n ,i,)neN d ; nes m the domain V{1-L ND ) and that 

2 d 



\*n,L(*)\ 2 < L „ 



We have 



HZ D ($ rh L,$ n ,L)= ! \V<S> n , L (x)\ 2 dx- [ | V$„, L (x) fdx. (4.35) 

Ja l Jn L 

As ($ n ,L)n6N d are the eigenfunctions of the Laplacian restricted to the box Al 
with the appropriate boundary condition [19], we get the following properties, 

• for any n, m G N d such that n ^ m, we have 

/ $u,l( x ) ■ ®m, L (x)dx = [ V$ n , L (x) ■ V$ m , L (x)dx = 0, (4.36) 

JA L JAl 

• for any n G N d , 

' \V<5> n , L (x)\ 2 dx = a n (L). (4.37) 



L 



Now using equation (I4.35P and taking into account (14.361) and (I4.37P we get 
that, 

\HZ D (® n>L ,® n , L ) -a n (L)\ < a n (L) ■ (j) d \Q L \ (4.38) 

and 

\nZ D (* n , L ,* m ,L)\ < {a n {L))l{a m {L))^ ■ £) d • |fi L |. (4.39) 



Applying Lemma 14.21 with 



a = E > 0, e 1 = {j) d - | Q L \< 1, e 2 = a • (j) d - \ Sl L \, 
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we get that there exist constants C\ , C 2 > such that 

N{H%f, G x ■ E) > C 2 ■ ■ L d . (4.40) 

The right hand side of equation (I4.40j) is due to the fact that for —A, there 
exists C 2 > such that [JS] 

%{n;a n (L)<E} = C 2 -E d / 2 L d . (4.41) 

The proof of Theorem 14.11 is ended by taking into account fl4.40p and the 
following equation [9] 

j- d .N(H» L D ,E)<N(E). (4.42) 

□ 

The proof of Lemma 14. 2b We start by recalling the following expression 
of the counting function given in terms of a supremum of the dimension of 
all linear subspace £ in the domain of A which we denote by T>. 

N(A, E) = sup{dim£\(p,A(p) < E(<p,<p), for all if G £}. (4.43) 
ecd 

When E\ < 1, (p 1} • ■ ■ , ip n span a subspace £ n of dimension n. For any $ G £ n 
there exist (non- unique) coefficients ci, • • ■ , c n G C such that 

n 

$ = ^gifi. 
i=l 

Using (I4.3ip . one gets 

ran n 

1=1 ij i=l 

and 

n n n 

(<M*><£ Y N ■ \ c k\ ■ {((PijAfy,) -anSijl < {a + e 2 )^2\ci\ 2 . 

i=l M=l i=l 

(4.45) 

OL -\- E 2 

The proof of Lemma [4.21 is ended by setting E = . 

1 - Si 
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4.2 Examples 

Let us consider some particular cases. 



4.2.1 The periodic case. 

Consider 

n= |J A j9 ( 7 ),/3<i, 



here A j g(7) is the cube of center 7 and side length /3. In this particular case we 
get \Ql\ = L d (3 d , and the assumption of Theorem 14.11 is satisfied for j3 < |. 

Remark 4.3. The periodic model includes random displacement models. 
4.2.2 The Anderson case. 

Consider now the case of the model described in subsection 12.11 which is the 
object of our studies. 

n= |J ^(5 + 7). 



In this situation we get: 

For P-almost all u and any a >| <S | -p and all L large enough we have 

7^l|<« (4.46) 

So the assumption of Theorem 14. II is satisfied for | S \ -p < ^. 
To get (14 .46 p . we use the fact that S is bounded, so there exists L such 
that S C Ai , and we get 



So 



= u 7eZ dw 7 (5 + 7) n A L c u 7eALo+L uj 7 (S + 7). 



p-l^l ^ J? £ <«Vl 5 l (4.47) 
7eAL 0+ L 

£ t'^i-'ct E «.)■ <*«) 
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Now by taking L big enough we get the result as the term between parenthisis 
in eqution (14.48 p . converges to E(wo) = V by the strong law of large numbers. 

4.2.3 The Poisson case. 

Let us starts by giving the model. Let m(dx) be a Poisson random measure 
on M. d . This means that if {A\,--- ,A n ) are pairwise disjoint Borel sets in 
M d , the random variables m(Ai), ■ ■ ■ m(A n ) are independent, and when A is a 
bounded Borel set, for k = 0, 1, • • • , the random variable has the distribution 



Here c is the the concentration defined by ~E(m(A)) = c ■ \ A\. If are the 
atoms of the Poisson measure then m can be written as 



¥{m(A) = k} = e 



MA\) k 



k\ 





In the above notation we have 



r w = ^ = |J(5 + 6). 



So we get that 




• m(A Lo +L)- 



(4.49) 



L d ' l ° l ' K (L + LY' 



As 




1 



(4.50) 
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we get that the assumption of Theorem 14.11 is satisfied for 

C-\S\ < -r. 
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